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Abstract 

We perform a phenomenological analysis of dijet production in 
double pomeron exchange at the Tevatron. We find that the CDF 
Run I results do not rule out the presence of an exclusive dijet 
component, as predicted by Khoze, Martin and Ryskin (KMR). 
With the high statistics CDF Run II data, we predict that an ex- 
clusive component at the level predicted by KMR may be visi- 
ble, although the observation will depend on accurate modelling 
of the inclusive double pomeron exchange process. We also 
compare to the predictions of the DPEMC Monte Carlo, which 
contains a non-perturbative model for the central exclusive pro- 
cess. We show that the perturbative model of KMR gives dif- 
ferent predictions for the di-jet Ex dependence in the high di-jet 
mass fraction region than non-perturbative models. 



1 Introduction 



The detection of new particles in central exclusive production at the LHC has received a 
great deal of attention recently (see for example [1-6] and references therein). By cen- 
tral exclusive production, we refer to the process pp ^ p ® 4> ® p, where © denotes 
the absence of hadronic activity ('gap') between the outgoing protons and the decay 
products of the central system 0, shown schematically in Fig. 1 (a). An example would 
be Standard Model Higgs boson production, for which the central system could consist 
of either 2 b-quark jets or the decay products of two W bosons, and no other activity. 
Khoze et al. calculate the cross section for the central exclusive production of a 120 
GeV Standard Model Higgs boson to be 3 fb at the LHC, falling to ~ 1 fb at M^^ = 200 
GeV [7, 8], and orders of magnitude larger for certain MSSM scenarios [9]. This is 
large enough to be observable, and there is currently a proposal to install forward proton 
detectors as a future upgrade to the LHC experiments [10]. It is therefore important to 
confront the KMR calculations with data as soon as possible, in particular by searching 
for high-rate central exclusive processes at the Tevatron. 

There are several exclusive processes which, according to the KMR calculations, should 
have large enough cross sections to be observable now. The two cleanest signatures 
are central exclusive xc meson production and central exclusive di-photon production. 
These processes are under study by both the CDF and D0 Collaborations, although at 
the time of writing, no results have been published. By far the highest rate process 
is predicted to be central exclusive di-jet production. Di-jet production in the inclu- 
sive process pp ^ p ® X ® p, where X can be any centrally produced system, was 
studied by the CDF Collaboration in the Tevatron Run I data [11]. This process, of- 
ten termed double pomeron exchange, is conventionally modelled as shown in figure 
1(b). In Regge-inspired language, it can be pictured as the emission of two pomerons 
(or reggeons if the longitudinal momentum loss of the protons is large) which collide to 
produce the central system X. Since the pomerons have a partonic structure, described 
by structure functions, there are necessarily always pomeron remnants; that is to say 
there are no truly exclusive events in this picture. 

Whilst from a theoretical perspective, the exclusive process is distinctly different from 
the inclusive, from an experimental viewpoint the definition of exclusive and inclusive 
is somewhat arbitrary. The approach adopted by the CDF collaboration is to define the 
quantity Rjj as the fraction of the invariant mass of the central system contained within 
the two highest Et jet cones. If the cone radius is large enough such that out of cone 
effects are small, then Rjj should be peaked at or near 1 for exclusive events since there 
will be no pomeron remnants. If the exclusive process is present with a large enough 
cross section relative to double pomeron exchange then an excess of events over that 
predicted by the Regge models should be observable at high Rjj. 

In this paper, we carry out a Monte Carlo simulation of di-jet production in double 
pomeron exchange, and compare directly to the Run I data. We use the POMWIG 
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Fig. 1: The exclusive production process (a) and central inclusive (or double pomeron) process (b). 

Monte Carlo generator to simulate the inclusive double pomeron exchange process [12], 
and ExHuME [13] to simulate the exclusive contribution, as calculated by KMR. An al- 
ternative model for the exclusive process, based on the calculations of Bialas and Land- 
shoff (BL) [14], is implemented in the DPEMC Monte Carlo [15]. There are two key 
differences between the BL and the KMR approaches. Firstly, in the KMR model, the 
two gluons at the proton vertex couple perturbatively to the off-diagonal unintegrated 
gluon distribution of the proton. In the BL approach, the proton vertex is treated as 
non-perturbative, and parameterised using a Regge-motivated ansatz. This leads to a 
significant difference in the rapidity distributions of the central system; in the BL ap- 
proach the rapidity distribution is much flatter, since the cross section has only a weak 
x-dependence, entering via the soft pomeron intercept, whilst in the KMR model there 
is a much stronger x-dependence from the gluon distributions in the protons. Secondly, 
the cross section falls much more quickly as the di-jet Et increases in the KMR model. 
This is a reflection of the fact that the phase space for gluon emission into the gap in- 
creases as Et increases, and therefore the exclusive cross section decreases. This effect 
is not present in the BL approach. As we shall see, this phase space effect has observable 
consequences which could allow the two models to be separated in the Tevatron Run II 
data. 

CDF have also recently presented high statistics preliminary results on diffractive di-jet 
production in Run II [16, 17]. Since the Run II results are only preliminary at the time of 
writing, we do not compare directly, but make a prediction based on the kinematic range 
quoted in [16]. Our aim is to assess whether the exclusive process should be visible in 
the CDF data. 



2 Monte Carlo simulation 

To simulate the exclusive process, we use both the ExHuME 1.3.1 generator [13], in- 
terfaced to Pythia 6.205 [19,20] for parton showering and hadronisation, and DPEMC 
2.4 [15]. To simulate the double pomeron exchange process, we use POMWIG 1.3 [12]. 
POMWIG simulates diffractive collisions using a Regge factorisation ansatz, with the 
parameters of the pomeron and reggeon flux terms and structure functions extracted 
from diffractive DIS data by the HI CoUaboration at HERA [21]. The POMWIG model 
generates only the inclusive process of figure 1(b) - i.e. there will always be pomeron 
remnants in the event. We use the default POMWIG parameters as described in [12], 
with one exception. The default POMWIG behaviour up to v 1 .2 was to treat the pomeron 
(and reggeon) remnants in the same way that HERWIG treats the photon remnants, i.e. 
the valence partons in the pomeron and reggeon are defined to be quark - antiquark pairs. 
POMWIG vl.3 introduces the option to define the valence partons in the pomeron to be 
gluons. It is this option that we take to be the default behaviour for the pomeron in 
this paper. For the Reggeon, we define the valence partons to be quarks, as in version 
1.2. We discuss the effects of the different remnant treatments in section 4. Pomeron 
and reggeon exchanges are generated separately, and added according to their gener- 
ated cross sections. POMWIG does not include a model of gap survival probability, 
and therefore over-estimates diffractive cross sections at the Tevatron and LHC. The ap- 
proach we adopt is to scale the final POMWIG prediction, after the experimental cuts 
and smearing described in section |3l to the CDF Run I measurement [11]. This as- 
sumes that the gap survival factor is independent of kinematic effects, which is probably 
a reasonable approximation given the level of accuracy we are working to here. In the 
exclusive case, the gap survival factor is a prediction of the KMR model. We use the 
ExHuME default value of = 0.045 at the Tevatron for both ExHuME and DPEMC 
^ The normalisation of the ExHuME cross section is particularly sensitive to the gluon 
distribution in the proton (which enters roughly speaking to the fourth power). We use 
the MRST2002 parton distribution functions [23] as the ExHuME default. To gain some 
estimate of the uncertainty, we also use the CTEQ6M pdfs [24] in our Run II predictions, 
which lead to a larger cross section prediction ^. 

3 Results 

3.1 CDF Run I results 

We first turn our attention to the published CDF Run I results [11]. CDF have proton 
tagging detectors only on the outgoing anti-proton side of the experiment. Tagging the 

'The authors of DPEMC suggest that the gap survival factor to be used in DPEMC at the Tevatron is 
0.1, although the uncertainty in the normalisation of the DPEMC cross section is at least a factor of 2, 
which covers the lower value we use here [22]. 

^For a recent analysis of the effect of using different parton distribution functions on the central exclu- 
sive cross section, see [25] 
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outgoing p restricts its longitudinal momentum loss to the range 0.035 < < 0.095, 
and the transverse momentum \tp\ < 1 GeV^. On the outgoing proton side, CDF have 
no proton taggers, so the analysis relies on the observation of a rapidity gap (i.e. no 
particles above a noise threshold) in the forward detectors, covering the range 3.2 < 
\r]\ < 5.9. This corresponds to an approximate range of proton longitudinal momentum 
loss 0.01 < < 0.03. Following the CDF analysis, jets are found using the cone 
algorithm with cone radius R — 0.7. In table 1 we show the generated cross sections 
from POM WIG, ExHuME and DPEMC in the CDF Run I kinematic range [11], defined 
as follows; 

-4.2 < r)jetsi,2 < 2.4 (2) 
0.035 < < 0.095 (3) 
\tp\ < IGeV^ (4) 
0.01 < < 0.03. (5) 

We note that, at the hadron level, POMWIG overshoots the measured cross section by a 

factor of ~ 8. This would imply an effective gap survival factor of 0.14. As we shall see 
below however, a combination of the smearing of the Monte Carlo, and the rapidity gap 
definition used by CDF to calculate the measured cross section, changes this result by a 
factor of ~ 2. 





ET,rmn (GcV) 


aroT (nb) 


(^Rjj>o.8 (nb) 


CDF 


7 


43.6±4.4±21.6 


< 3.7 


POMWIG 


7 


320.74 


0.69 


ExHuME 


7 


0.48 


0.08 


DPEMC 


7 


1.32 


0.26 


CDF 


10 


3.4±1.0±2.0 


n/a 


POMWIG 


10 


49.74 


0.36 


ExHuME 


10 


0.25 


0.07 


DPEMC 


10 


0.68 


0.22 



Table 1 : Generator Level Cross Sections and CDF Data for Run 1 . The POMWIG result has no effec- 
tive survival factor applied as this is calculated using the smeared data. Rjj is calculated by using the 
generated ^1,2 values and unsmeared final state hadrons. 

The CDF Collaboration results are not fully corrected for detector effects. In order to 
compare directly to the published CDF distributions we must therefore smear the energy 
and momenta of the final state particles from the Monte Carlos. This is done according 
to the figures quoted in the CDF technical design reports [26, 27] before jet finding. 

The first step in the CDF analysis is to define a single diffractive di-jet data set, using the 
proton tagging detector on the outgoing anti-proton side, in the kinematic range defined 
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by cuts (1), (2), (3) and (4) above. We smear the energy and momenta of the final state 
particles, but take the and tp values from the Monte Carlo unsmeared. A calorimeter 
noise suppression cut is applied at Eceii < 100 MeV. We apply this cut to the smeared 
energies of the final state particles. 

A rapidity gap on the outgoing proton side is required by demanding no particle hits 

in the Beam-Beam Counters (BBC), covering the pseudorapidity region 3.2 < rj < 5.9, 
and no hits in the Forward Calorimeter towers (FCAL), 2.4 < \ri\ < 4.2, with E'^e// > 1-5 
GeV. This cut approximates to the range given in (5). We apply the above rj cuts and 
calorimeter noise thresholds to the smeared final state Monte Carlo particles. The value 
of is then calculated from the calorimeter information by 



Cp=^ H Ere^ (6) 

V particles 

where the sum is over all final state particles excluding the outgoing protons. The 

calorimeter-based measurement is calibrated by comparing the direct ^p measurement 
from the proton tagging detector on the outgoing anti-proton side with the corresponding 
calorimeter-based measurement of ^p, 

Cp—^ E Ere-^ (7) 

V particles 

This method allows CDF to calculate an average correction factor which can be used to 
correct the calorimeter-based ^p measurement. CDF find that the ^p value measured in 

the calorimeter should be multiplied by a correction factor of 1.7. We have performed 
an identical procedure on the smeared Monte Carlo and found a smaller correction (1.2) 
is required. 

This could be due to the fact that CDF had very low statistics in the data sample from 
which the correction factor was derived, or could indicate that extra detector effects, 
apart from smearing, are present. It could also of course be due to inaccuracies in the 
POMWIG model. We use our correction factor of 1.2 in the remainder of the Run I 
analysis - i.e. we make a correction for the kinematical reconstruction method of ^p, but 
allow for the possibility that additional effects might have been corrected for in the CDF 
data. 

In figure 2 we compare the ^p and ^p distributions from POMWIG and ExHuME to 
the CDF Run I data. The agreement is reasonable in ^p, which is measured directly in 
the proton tagger. For the calorimeter-based ^p distribution, we reproduce the overall 
shape of the data. We attribute the differences in certain bins to statistical fluctuations in 
the data, or to detector effects that we are unable to reproduce, since one would expect 
this distribution to be smooth, except in the lowest { region, where kinematic factors 
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cause the distribution to turn over. In figure 3 we show the mean jet Et and mean 
jet distributions of the two highest Et jets compared to the CDF Run 1 data. We 
obtain reasonable agreement. We attribute the undershooting of the data in the lowest 
Et bin, and the more peaked mean rj distribution in the Monte Carlo, to detector effects. 
The DPEMC samples (not shown) give similar results. The distributions are reasonably 
well described by POMWIG alone. This was also the conclusion of a similar analysis 
presented in [28]. 



p n 111,1 

0.04 0.05 0.06 0.07 0.08 0.09 



- POMWIG + ExHuME 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.0! 



Fig. 2: The and distributions compared to CDF Run I data [11]. is measured directly in the 
forward proton tagging detector, and ^p is measured using the central detectors, as described in the text. 
The solid line the sum of the ExHuME and POMWIG predictions. 

To calculate the cross section, CDF restrict the DPE sample to the range 0.01 < < 
0.03. The double pomeron exchange cross section is obtained by taking the ratio of the 
number of double diffractive di-jet to single diffractive di-jet events, and multiplying 
by the single diffractive di-jet cross section. Using this method, CDF measure the DPE 
di-jet cross section, in the kinematic range defined by cuts (1) - (5), to be a^^^ — 
43.6 ± 4:A{stat) ± 21.6{syst) nb. 

With the DPE cross section defined in the above kinematic range (i.e. cutting on the 
generated values of ^p, the scaled calorimeter measurement of ^p, and the smeared final 
state particles in the jet finding), we find that POMWIG over-estimates the measured 
cross section by a factor of ~ 4, leading to an effective gap survival factor of 0.27. 
This fixes the normalisation of POMWIG. The discrepancy of a factor of ~ 2 between 
the hadron and detector level normalisation factors occurs because of the rapidity gap 
requirement, which results in the loss of events in the 0.01 < < 0.03 range and not 
just at high ^p. We note that POMWIG does not simulate reggeon - pomeron collisions. 
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and may therefore underestimate the cross section at the Tevatron as predicted from the 
HERA data. The gap survival factor should therefore be considered as an 'effective' gap 
survival factor which takes into account such effects. It is also worth noting that there is 
a systematic error of ~ 50% on the overall normalisation of the CDF cross section. 




Fig. 3: The mean jet and mean jet 77 distributions of the 2 highest Ex jets, defined as described in the 
text. The data points are the CDF Run I resuhs [11]. The dashed line is the prediction of the ExHuME 
Monte Carlo, and the solid line the sum of the ExHuME and POM WIG predictions. 

The di-jet mass fraction Rjj is defined as Rjj = Mjj/Mx, where Mx is the invariant 
mass of the diffractive central system and Mjj is the invariant mass of the 2 highest Et 
jets. Experimentally, the total mass in the calorimeter is a difficult quantity to measure. 
CDF reconstruct Rjj as follows (we follow this procedure, using the smeared final state 
particles from the Monte Carlos) 



MIj = X1X2S (8) 
^1 = ^ E ^Te" (9) 

V "5 jetsl,2 

X2 = ^ ^Te-^ (10) 

V jetsl,2 

where the Xi are approximately, in the POMWIG model, the fraction of the pomeron 
momenta carried by the struck partons in the pomerons (often termed /?). Similarly, 

Mi = ei6s (11) 
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Fig. 4: The fraction of the invariant mass of the central system contained within the 2 highest Et jet 
cones. The dashed Une is the prediction from ExHuME (i.e. the central exclusive process), the dotted line 
is the prediction from POMWIG, and the solid line is the sum of the ExHuME and POMWIG predictions. 
The output of both Monte Carlo generators was smeared to simulate detector effects, as described in the 
text. The data points are the CDF Run I results [11]. 

ep = ^ E ^^e" (12) 

V particles 

= ^RP (13) 

where the sum is over all final state particles (excluding the outgoing protons), and ^rp 
is the fractional longitudinal momentum loss of the outgoing anti-proton measured in 
the Roman Pot. We apply the correction factor of 1.2 to as discussed above. 

In figure 4 we show the smeared Rjj distributions compared to the CDF Run I results. 
Between Rjj = 0.2 and Rjj = 0.7, POMWIG alone gives a reasonable description of 
the data. This is consistent with the findings of [28]. The POMWIG prediction is shifted 
to higher Rjj values however. There are two related reasons for this shift. Increasing the 
correction factor applied to results in a shift to lower Rjj (since doing so increases the 
measured Mx value in the denominator of Rjj). As we noted above, the CDF correction 
factor of 1.7 had a large error due to limited statistics, and may therefore have been too 
large. Similarly, the distribution shown in figure 2 rises more steeply with increasing 
in the data than in POMWIG, although the statistical errors on the data are large. This 
again has the effect of increasing the Mx measurement in the data, and shifting the Rjj 
distribution to lower values. The approach we take here is to attribute these shifts to 
statistical fluctuations and detector effects, which we are unable to simulate. 

For Rjj > 0.8, POMWIG generates few events, and the distribution becomes domi- 
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nated by exclusive events generated by ExHuME (or DPEMC). CDF had insufficient 
statistics to probe this region in Run I. In table 2 we show the cross sections predicted by 
POMWIG, ExHuME and DPEMC after detector smearing, in the range 0.01 < <0.03 
using the effective gap survival of 0.27 for POMWIG, and the predicted gap survival 
factor of 52 = 0.045 for ExHuME and DPEMC. 





Ermin (GeV) 


f^TOT (nb) 


7>o 8 (nb) 


CDF 


1 


43.6±4.4±21.6 


< 3.7 


POMWIG 


7 


42.53 


0.01 


ExHuME 


7 


0.59 


0.03 


DPEMC 


7 


1.29 


0.09 


POMWIG + ExHuME 


7 


43.12 


0.04 


POMWIG + DPEMC 


7 


43.8 


0.10 


CDF 


10 


3.4±1.0±2.0 


n/a 


POMWIG 


10 


6.91 


<0.01 


ExHuME 


10 


0.28 


0.03 


DPEMC 


10 


0.60 


0.09 


POMWIG + ExHuME 


10 


7.19 


0.03 


POMWIG + DPEMC 


10 


7.52 


0.09 



Table 2: The cross section predictions from POMWIG (with an effective gap survival factor S*^ = 0.27), 
ExHuME and DPEMC, in the kinematic range described in the text, with detector smearing included. 
Also shown are the CDF Run I published cross sections, taken from [11]. 

DPEMC predicts a larger cross section than ExHuME after the CDF cuts. This is partly 
due to the soft ^ dependence of DPEMC which leads to more events being generated 
in the large { region. It should always be remembered however, that there are large 
uncertainties in the normalisations of the ExHuME and DPEMC predictions, and both 
the absolute and relative normalisations can vary within at least a factor of 2. Neither 
model is ruled out by Run I results. 



3.2 CDF Run II predictions 

Having used the published Run I data to fix the normalisation of the POMWIG sample 
(5^ = 0.27), we can predict Rjj for the increased statistics (and differing kinematic 
range) of Run II, following the approach used in the CDF preliminary central di-jet 
production analysis described in [16]. The kinematic range for Run II is as follows; 

> 10 GcV (14) 

IVjetslM < 2.5 (15) 

0.03 < < 0.1 (16) 

< 0.1 (17) 
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tp\ < IGeVl 



(18) 



The larger range is due to a new rapidity gap definition, partially imposed by the new 
configuration of forward detectors at CDF Run 11 [18]. In table 3 we show the hadron 
level cross section predictions from POMWIG, ExHuME and DPEMC. We also show 
the ExHuME prediction using the CTEQ6M parton distribution functions, which lead to 
a larger cross section (although still significantly lower than the DPEMC prediction). 

In figure 5 we show the hadron level Rjj distribution predicted by POMWIG, ExHuME 
and DPEMC. Due to the effects of parton showering and hadronisation, neither model 
predicts a visible excess of events at large Rj j. In fact the majority of exclusive events 
lie below Rjj = 0.8. The DPEMC exclusive cross section is larger than that of Ex- 
HuME, as in the Run I kinematic range, although the caveats about the absolute and 
relative normalisations of ExHuME and DPEMC still apply. 

In figure 6 we show the Et distribution of the second highest Et jet in the 'exclusive' 
region Rjj > 0.8, for ExHuME and DPEMC. Here, there is a difference between the 
predictions of the models. As discussed in section 1, this is because the ExHuME (KMR) 
cross section falls rapidly as the di-jet Et rises, due to the increased phase space for 
gluon emission. In non-perturbative models, this suppression is not present, and the jet 
Et dependence of the cross section is therefore much flatter. In figure 6 we also show 
the jet Et dependence of events generated by ExHuME, DPEMC and POMWIG alone, 
normalised such that all curves pass through the same point at Et = 10 GeV, making 
the difference in the predicted slopes easier to see. In common with DPEMC, there is 
less suppression of the production cross section with increasing Et in POMWIG. Note 
also that the parton distribution functions used in ExHuME affect the Et distributions 
in a non-negligible way. CTEQ6M produces a softer Et dependence than MRST2002, 
although both are significantly steeper than the non-perturbative models, for the reasons 
discussed above. 

4 Summary and Discussion 

We have shown that the published CDF Run 1 double pomeron exchange data do not 
rule out the presence of an exclusive component at the level predicted by Khoze, Martin 
and Ryskin (KMR). The DPEMC prediction is also not ruled out by the Run 1 results. 
Looking forward to the Tevatron Run II data, we have made a prediction for the di-jet 
mass fraction, Rjj, based on the POMWIG Monte Carlo, which we scale to the Run 1 
measured cross section, and the absolute prediction of the KMR model from ExHuME. 
At the time of writing, only preliminary results are available for the high - statistics Run 
11 data. In the absence of final, published results, we present hadron level predictions, 
in the kinematic range used in the preliminary analysis [16]. As is clear from figure 

5 and table 3, POMWIG generates fewer events at Rjj > 0.8 than ExHuME. With 
sufficient statistics, and if detector effects are corrected for, we might therefore expect 
the Rjj distribution in the data to be higher than the POMWIG prediction at Rjj > 0.8. 
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(Ttot (.no; 


(^Rjj>o.8 (nb; 


rUMWlLi 


1 n 
iU 


1 OO 1 

ioo.iO 


U.iU 


E/Xriuiviij/ t^ivirvo 1 zyjyjz ) 


1 n 


U.oZ 


U.zo 








r\ /1 7 
U.4j 




1 n 


9 61 
Z.Di 


U.oU 


r UiVl W LKj + liXrlUiVlrl, (^iVlKo i zUUZ ) 


1 r\ 
iU 


1 QQ QQ 




I'UMWlLr + ilXrluMil (L- 1 ilV^oM) 


iU 


ioy.oi 


yj.jj 




iU 


lyu. / / 


u.yu 


rUMWlLr 


oc 

ZJ 


A C\Af\ 


O.vOo 




Z J 


U.Uivj 


U.UiZ 


ExHuME (CTEQ6M) 


25 


0.037 


0.027 


DPEMC 


25 


0.176 


0.118 


POMWIG + ExHuME (MRST2002) 


25 


0.956 


0.020 


POMWIG + ExHuME (CTEQ6M) 


25 


0.977 


0.035 


POMWIG + DPEMC 


25 


1.116 


0.126 



Table 3: The cross section predictions from POMWIG (with an effective gap survival factor S"^ = 0.27), 
ExHuME and DPEMC, in the CDF Run II preliminary kinematic range as described in the text. Also 
shown is the ExHuME prediction using the CTEQ6M proton pdfs. 

This could be interpreted as a sign of an exclusive component in the data, given that 
POMWIG is able to adequately model the high Rjj region. This is an assumption, 

but we believe it to be a reasonable one, for the following reason. In order to produce 
more events at large Rjj, the pomeron structure function would have to be larger at 
high-/? than it is in POMWIG. In figure 7, we compare the effect of varying the structure 
function by comparing the POMWIG default (HI fit 2 [21]) to the HI fit 3 peaked gluon. 
Even with the relatively extreme high-/? behaviour of the HI fit 3 gluon distribution 
(which is certainly too large), an excess at high Rjj would still be present in the data 
above the POMWIG prediction if an exclusive component at the level predicted by KMR 
is present, since the cross section from ExHuME alone is as large or larger than the 
POMWIG HI fit 3 prediction in the highest Rjj bins. As discussed above, there is an 
uncertainty in the absolute magnitude of the KMR prediction, which is largely due to 
the uncertainty in the knowledge of the gluon distributions of the proton. In figure 7 
we also show the ExHuME prediction using the CTEQ6M parton distributions (see also 
table 3). In this case, the excess over the POMWIG prediction is larger. 

The treatment of the pomeron remnant can also have an effect on the POMWIG predic- 
tions at high Rjj. As discussed in section 2, POMWIG vl.3 introduces the option of 
defining the pomeron valence partons to be either quarks or gluons. Since in the HI pic- 
ture, the pomeron is predominantly a gluonic object, the more natural choice is to define 
the valence partons to be gluons. In figure 8 we compare the Rjj distributions when 
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Fig. 5: The Rjj distributions at the hadron level predicted by POMWIG + ExHuME (left hand plot) and 
POMWIG + DPEMC (right hand plot), in the CDF Run II kinematic range, as described in the text. 

the pomeron valence partons are defined to be quarks (the POMWIG default behaviour 
prior to vl.3) or gluons. The POMWIG prediction is significantly lower at high Rjj for 
the quark case. This is because in di-jet production, the parton entering the hard scatter 
from the pomeron is more likely to be a gluon. If the valence partons in the Pomeron 
are defined to be quarks, HERWIG will be forced to continue parton showering until 
a quark-antiquark pair is produced. This results in increased activity in the remnant 
region, and hence high Rj j values are less likely. 

Both the DPEMC and POMWIG models have significantly different jet Et dependences 
to ExHuME (with both the CTEQ6M and MRST2002 partons distributions), as shown 
in figure 6. The steeper jet Et dependence seen in ExHuME is a direct result of the fact 
that the KMR model takes account of the increasing phase space for gluon emission as 
Et increases. If, on the other hand, the high Rjj events are due to the tail of a POMWIG 
(or a DPEMC) - like model, then a flatter Et dependence should be observed. The jet 
Et distribution in the high Rjj region is therefore a possible test for the presence of a 
perturbative exclusive component (i.e. KMR) in the double pomeron di-jet data. 
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Fig. 6: The Et distribution of the second highest Et jet in the region Rjj > 0.8. The predictions of 
POMWIG + ExHuME (MRST2002), POMWIG + ExHuME (CTEQ6M) and POMWIG + DPEMC are 
shown in the left hand plot. The kinematic range is as described in the text. In the right hand plot, the 
predictions of ExHuME alone (with MRST2002 and CTEQ6M structure functions), POMWIG alone and 
DPEMC alone are shown with the curves normalised such that they all pass through the same point at 
Et = 10 GeV. 
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